Field dependence of the vortex core size in a multi-band superconductor 
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The magnetic field dependence of the vortex core size in the multi-band superconductor NbSe2 
has been determined from muon spin rotation measurements. The spatially extended nature of the 
quasiparticle core states associated with the smaller gap leads to a rapid field-induced shrinkage of 
the core size at low fields, while the more tightly bound nature of the states associated with the 
larger gap leads to a field-independent core size for fields greater than 4 kOe. A simple model is 
proposed for the density of delocalized core states that establishes a direct relationship between the 
field-induced reduction of the vortex core size and the corresponding enhancement of the electronic 
thermal conductivity. We show that this model accurately describes both NbSe2 and the single-band 
superconductor V3SL 

PACS numbers: 74.25.Ha, 74.25.Jb, 74.25. Qt, 74.70.Ad, 76.75.-K 



The discovery of MgB 2 Q, a seemingly conventional 
type-II superconductor with a high transition tempera- 
ture T c = 39 K, has sparked considerable experimental 
and theoretical interest in the physics of multi-band su- 
perconductivity (MBSC). MBSC, which was first treated 
theoretically in 1959 using BCS theory 0], implies that 
below T c distinct energy gaps open up on different sheets 
of the Fermi surface. In MgB2, a second smaller gap is 
induced on the ir bands through coupling to the intrinsi- 
cally superconducting a bands. While the high value of 
T c and the large boron isotope effect in MgB 2 are now be- 
lieved to originate from strong electron-phonon coupling 
on the two-dimensional a bands, it is also clear that sev- 
eral properties of MgB2 are a manifestation of MBSC 
EJ El IE El ■ Large vortex cores have been observed on 
the 7r bands by scanning tunneling spectroscopy [3J ■ Mi- 
croscopic calculations of the local density of states (DOS) 
in a two band superconductor predict that the smaller 
gap will give rise to spatially extended quasiparticle (QP) 
states at low field [a 13- At higher field, these loosely 
bound states become delocalized, with the core size now 
being determined by the more localized states associated 
with the larger gap. This model is supported by the mag- 
netic field dependence of the electronic thermal conduc- 
tivity 0,0] and the specific heat 0,0,0], techniques 
which are sensitive to the delocalization of QPs. 

Muon spin rotation (yitSR) has proven to be an ac- 
curate method for determining the field dependence of 
the vortex core size [15J, however single crystals are re- 
quired for this kind of measurement. Although there have 
been several studies of MgB 2 in the vortex state by /iSR 
00000, these have bee n done on powdered 
samples because of the difficulty of growing large, clean, 
MgB 2 single crystals. 

NbSe 2 has long been thought of as a typical s-wave 
type-II superconductor. However, there is now convinc- 
ing evidence that it too is a multi-band superconductor 
pll I22I 0, 0| . Experiments suggest that an energy gap 
of magnitude ~ 1 meV exists on two bands derived from 
Nb Ad orbitals, and a gap ~ 2-3 times smaller also exists 



on a band derived from Se Ap orbitals. In the quest to 
fully understand the physics of MBSC, it is important 
to study superconductors other than MgB 2 , and to this 
end we performed fiSK measurements of the field depen- 
dence of the vortex core size in NbSe 2 at very low tem- 
perature. A distinct advantage of NbSe 2 is the availabil- 
ity of large, high quality, single crystals. Single crystals 
allow for the alignment of the initial muon spin polar- 
ization along a single principal crystallographic axis. In 
this case the /xSR lineshape resembles the asymmetric 
theoretical internal magnetic field distribution n(B) of 
the vortex lattice. In particular, the high-field tail is the 
contribution of the vortex cores to n(B) (see Fig. 1) 0|. 
Previous /xSR measurements of the field-dependence of 
the core size in NbSe 2 were performed at T>2.3 K [26| . 
At such temperatures, thermal excitations of the bound 
core states (the so-called 'Kramer-Pesch effect' 0,0) 
result in significant overlap of the QP wavefunctions of 
adjacent vortices. When this happens the QPs delocal- 
ize, giving rise to a strongly field-dependent core size [29j, 
as predicted by the microscopic theory |30| . 

In order to freeze out thermal excitations of the bound 
core states and fully isolate the effect of magnetic field 
on the core size, we have carried out /iSR measurements 
on a single crystal of NbSe 2 (T c = 7.0 K and upper 
critical field H C 2 = 45 kOe) in a dilution refrigerator 
at T = 20 mK. This also permits a direct comparison 
with electronic thermal conductivity K e data, which is a 
measure of the extent of QP delocalization. The experi- 
ment was performed at the Tri-University Meson Facility 
(TRIUMF), Vancouver, Canada, with the magnetic field 
applied parallel to the crystallographic [001] axis. For 
each measurement, the field was applied at T > T c be- 
fore cooling the sample to the desired temperature. 

The spin-polarized positive muons implanted into the 
sample stop randomly on the length scale of the vor- 
tex lattice, and therefore uniformly sample n(B). Each 
muon precesses around the local field B at its site at the 
Larmor frequency lu — j^B, where 7 M /2-7r = 135.5342 
MHz/T is the muon gyromagnetic ratio. On decay of 



2 




54 55 56 

Frequency (MHz) 



57 



FIG. 1: Fast Fourier transforms of both the time dependence 
of the muon spin polarization in NbSe2 at T = 20 mK and 
H = 4 kOe (open circles) and the fit (solid curve) to Eq. (1). 
Inset: muon spin precession signal and fit to Eq. (1) (in a 
reference frame rotating at 53.3 MHz). 



the muon after an average lifetime of 2.2 fis, a positron 
is emitted preferentially along the direction of the muon 
spin. The time evolution of the muon spin polarization 
is determined by detecting decay positrons from an en- 
semble of ~ 2 x 10 7 muons. The exact functional form 
of the muon spin polarization depends on n(B). Further 
details of the exp erimental technique used here can be 
found in Ref . |25| . The fiSR time spectra were fit assum- 
ing a Ginzburg-Landau (GL) model for the spatial field 
profile given by 



B(r)=B (l-6 4 )]T. 



u Ki(u) 



G 



(1) 



where b = B/B C 2, Bq is the average internal field, G are 
the reciprocal lattice vectors, K\(u) is a modified Bessel 
function, u 2 = 2^ b G 2 (l + b 4 )[l - 26(1 - 6) 2 ], £ a6 is the 
GL coherence length, and X a b is the magnetic penetra- 
tion depth. As explained in Ref. |15| , previous fiSK works 
have demonstrated that the field-dependence of the pa- 
rameter £ a ;, reflects changes in the vortex core size due 
to changes in the electronic structure of the vortex cores. 

Figure 1 shows fast Fourier transforms (FFTs) of both 
the muon spin precession signal at H = 4 kOe and the 
fit to Eq. p. The FFTs closely resemble n(B) for a 
hexagonal vortex lattice, but are significantly broadened 
by the apodization used to reduce ringing and noise in 
the FFT due to the finite time-range and the short muon 
lifetime, respectively. The lineshape is also broadened by 
a Gaussian distribution of fields from randomly oriented 
nuclear dipolar moments. The width of this distribution 
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FIG. 2: Field dependence of the vortex core size (, a b at T = 
20 mK (squares) and the electronic thermal conductivity k s 
normalised to its normal-state value tt n (circles). The dashed 
lines are guides for the eye. 



was determined to be « 0.25 MHz. In addition, there is a 
background peak at ~ 54.1 MHz originating from muons 
stopping outside the sample. 

Recently, /iSR measurements on powder samples of 
Mgi_ a; Al :z 32 were modeled assuming two distinct energy 
gaps, and a two-component field distribution for the vor- 
tex lattice [l^. Specifically, this model assumes two dis- 
tinct coherence lengths, which implies that both large 
and small vortex cores exist simultaneously in the sam- 
ple. However, for NbSe2 we find that the data are well 
described by Eq. (JJJ . A visual inspection of Fig. 1 shows 
that the fit captures all of the main features of the /xSR 
lineshape. This is consistent with recent theoretical work 
predicting that the formation of two vortex sublattices is 
energetically forbidden in a two-gap superconductor |3l| . 

Recent ^SR measurements of the field dependence of 
the core size in the single- band superconductor VaSi |29| 
have confirmed theoretical predictions of field-induced 
core shrinkage due to QP derealization [3(| • In the case 
of NbSe2, the derealization of core states is dependent 
on both energy gaps. At low fields the rapid reduction 
of Cafc and the simultaneous rapid increase of n e with in- 
creasing magnetic field are consistent with this regime 
being dominated by loosely bound states associated with 
the smaller energy gap H, . If the loosely bound core 
states associated with the smaller gap have completely 
delocalized by H w 4 kOe, then we expect the behav- 
ior of at larger fields to be dominated by the larger 
gap. The saturation of above 4 kOe indicates that 
the QPs are highly localized in the vortex cores at higher 
fields. This is in line with STM measurements which 
show that the bound core states at H = 10 kOe are more 
spatially confined than those at H = 1 kOe j^J ■ We note 
that even though the shrinking of the core size saturates 
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FIG. 3: Field dependence of the magnetic penetration depth 
in NbSe 2 at T = 20 mK. 



above 4 kOe, K e still exhibits a considerable dependence 
on field. The reason is that although the number of de- 
localized QPs per vortex is no longer changing, the in- 
creasing number of vortices increases the total number of 
delocalized QPs (i.e. heat carriers) in the sample. 

Given that the core size at low (high) field is primarily 
determined by the small (large) gap, one might naively 
expect that the BCS relationship £ ~ v f /tt/S.q 1 where v F 
is the Fermi velocity and Ao is the magnititude of the 
energy gap, could be used to determine the ratio of the 
high and low field core sizes. However, recent calculations 
based on GL theory have shown that this relationship is 
not applicable to a multi-band superconductor and that a 
simple expression for the ratio of the two core sizes only 
exists in the limit of zero interband coupling [33^. For 
finite interband coupling, the ratio must be calculated 
numerically and is dependent on a number of material 
parameters. 

As can be seen from Eq. (1), the magnetic penetration 
depth X a b is also a fitting parameter in our analysis, ex- 
tracted simultaneously with Fig. 3 shows that X a b 
increases strongly with H at low fields and displays a 
weaker, though significant, dependence at higher fields. 
The Volovik effect |34| has previously been invoked to ex- 
plain the field dependence of A ao in s-wave type-II super- 
conductors [3{| . This effect involves a shift of the quasi- 
particle energy spectrum by an energy ~ vp ■ vs , where 
vs is the supercurrent velocity, and it plays a signifi- 
cant role in superconductors with gap nodes. However, 
for an isotropic s-wave energy gap, a significant depen- 
dence of A a f, on H is only expected if the thermal energy 
ksT (where ks is Boltzmann's constant) is comparable 
to the magnitude of the energy gap. At T — 20 mK, 
ksT sa 1CP 3 meV, whereas the size of the small energy 
gap in NbSe2 is « 10 _1 meV. Therefore we cannot at- 
tribute the observed behaviour of X(H) to the Volovik 



effect. 

Calculations of the electronic structure of a single vor- 
tex using the Bogoliubov-deGennes equations have shown 
that the presence of extended quasiparticle states sig- 
nificantly modifies the current density, and hence B(r), 
around the vortex core [3g. With increasing magnetic 
field these modified regions overlap. Since Eq. (1) does 
not account for this, the effect on B(r) shows up in our 
measurements as a field dependent X a b- In other words 
X a b in Eq. (1) is an effective magnetic penetration depth. 
We note that the field dependence of X a b is stronger at 
low field due to the more rapid delocalization of QPs from 
the vortex cores. 

Next we propose a simple model which directly relates 
our measurements of the vortex core size to the measure- 
ments of the electronic thermal conductivity. In the con- 
ventional picture of a type-II superconductor, the density 
of localized QP states at the Fermi energy N(E F )\ oc is 
proportional to n£, 2 H, where 7r£ 2 is the area of a single 
vortex core and the density of vortices in the sample is 
proportional to H (37], |3^| . If we assume that the field- 
induced density of delocalized QP states is equal to the 
reduction in N(E F )\ OCl then 



N{E F ) dcloc <x (w$ - < 2 )ff, 



(2) 



where £o ~ 97 A is the low-field value of the core size. 
Note that (nQ — ir£, 2 )H is also proportional to the re- 
duction of the total core area in the sample due to the 
delocalization of QP core states. 

In a metal, n e = (l/3)Cv F l = (l/9)n 2 N(E F )v F lk 2 B T, 
where C is the heat capacity per unit volume, fee is Boltz- 
mann's constant, v F is the Fermi velocity, and I is the 
electron mean free path [39| ■ Since it is the delocalized 
QPs that carry heat then n e oc N(E F )deioc, and it follows 
from Eq. (2) that 



(3) 



To show that this relationship is physically valid, 
in Fig. 4 we plot the measured quantities n e and 
«o - n i 2 ) H against (H c2 {0)/ H) 1 ' 2 , which is pro- 
prtional to the intervortex spacing. We also plot our 
earlier measurements on VaSi (Ref. 29] ) in this way. It 
can be seen that at large intervortex spacing (low H) nei- 
ther quantity changes significantly. However, as the field 
is increased and the vortices are brought closer together, 
both quantities exhibit the same field dependence within 
experimental uncertainty. This lends further support to 
the conclusion that in both single-band VaSi and multi- 
band NbSe2 the field dependences of £ a t, and n e have a 
common underlying cause, namely the delocalization of 
bound QPs. 

In summary, we have measured the field dependence 
of the vortex core size in the multi-band superconductor 
NbSe2 at T = 20 mK. The observed field dependence is 
explained by the effects of two superconducting energy 
gaps on the bound core states. In addition, we have 
experimentally established a direct correlation between 
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the field-induced reduction of the vortex core area and 
the electronic thermal conductivity in both single-band 
and multi-band superconductors. 



FIG. 4: Total reduction in core area (left axis, solid squares) 
and electronic thermal conductivity K e (right axis, open tri- 
angles) plotted against (Heil H) 1 ^ 2 (which is proportional to 
the intervortex spacing) for V3S1 and NbSe2. Both quantities 
are normalised to their values at H ~ 70 kOe for VaSi and 
H ~ 10 kOe for NbSe2. The dashed lines connect the K e data 
points. (The VaSi data are shifted upwards by 0.5 on the 
vertical scale for clarity.) 
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